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ABSTRACT 



Given a constant of motion for the one- dimensional harmonic oscillator 
with linear dissipation in the velocity, the problem to get the Hamiltonian 
for this system is pointed out, and the quantization up to second order in the 
perturbation approach is used to determine the modification on the eigenvalues 
when dissipation is taken into consideration. This quantization is realized 
using the constant of motion instead of the Hamiltonian. 
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1. Introduction 



To find the Lagrangian and the Hamiltonian from the equations of motion of 
a given system (so called "inverse problem of the calculus of variations" or 
"the inverse problem of the mechanics") is not so trivial as one could think 
at first sight, even for one-dimensional systems and despite their existence is 
guaranteed [1] here. If the system is autonomous (the forces do not depend 
explicitly on time) or nonautonomous (total force depend explicitly on time), 
it is natural to try to find Lagrangians and Hamiltonians which do no depend 
explicitly on time or which depend explicitly on time for the latter. One pos- 
sible approach to find these quantities for a one-dimensional system is to find 
first a constant of motion of the system, and then to obtain the Lagrangian 
and Hamiltonian [2] . This constant of motion, of course, is chosen to be inde- 
pendent explicitly on time for autonomous systems, and explicitly depending 
on time, for nonautonomous systems. It has been shown that when one tries 
to get quantities (constant of motion, Lagrangian or Hamiltonian) which are 
explicitly depending on time for autonomous systems, there may be a concern 
about ambiguities [3] or about consistence on the mathematical quantities [4]. 

However, one common procedure that has been done is to guess an explicitly 
time depending Hamiltonian for the harmonic oscillator with exponential time 
dependence of its mass [5] which brings about the one-dimensional harmonic 
oscillator with linear dissipation in the velocity. Indeed, there has been a lot 
of studies on the quantization of this Hamiltonian [6] and decoherence of this 
quantum system [7]. However, since the one-dimensional harmonic oscillator 
with linear dissipation on the velocity is an autonomous system, one should 
get time independent dynamical quantities (constant of motion, Lagrangian 
or Hamiltonian) to describe it. After this, one should proceed to make the 
quantization of this system in a consistent way. In fact, an explicitly time 
independent constant of motion for this system has been found already [8]. So, 
in this paper we will study the classical damping behavior of the system with 
this constant of motion, and we point out the difficulties to get its associated 
Hamiltonian. Finally, we study its quantization using the same constant of 
motion and the velocity operator. This last study is done up to second order 
in perturbation theory and for weak dissipation. 
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2. Constant of motion and Lagrangian 



The one- dimensional harmonic oscillator with linear dissipation on the velocity 
is described by the equation 

mx + ax + kx = , (1) 

where m is the mass of the particle, a is the dissipative constant, k is the 
spring constant, x is the particle position, and x and x are its first and second 
differentiation with respect the time. By defining the new variable x — v, Eq. 
(1) can be written as the following autonomous dynamical system 



x = v 



and 



v = —u) 2 x — 2uj a v , 



(2) 



where u — yk/m is the natural angular frequency of the harmonic oscillator 
without dissipation, and u a = a /2m is the new dissipation parameter. A 
constant of motion for the system (2) is a function K = K(x, v) which satisfies 
the following partial differential equation [9] 



dK 

dx 



— (u 2 x + 2u a v) 



dK 

dv 







(3) 



The solution, K a , of this equation which satisfies the following limit 
lim Q ^ K a = mv 2 /2 + muj a x 2 /2 is given by [8] 



K a (x, v) = — I v + 2u a xv + u x J e 



2a; a G(i)/a;,ti;,ti) a ) 



(4a) 



where the function G has been defined as 



G(v/x,u,u a ) = < 



1 



2 Jul - u 2 



ui a + v/x 

1 



In 



u a + v/x - Jul - UJ 



u a + v/x + ^Ju 2 a - J 1 



arctan 



LU a + v/x 



to 2 - tul 



if uj 2 < uo 2 



iiu 2 = ul (45) 



if uj 2 > uj 2 



For weak dissipation, one has the following expression for the constant of 
motion 



K 



1 2 1 22 muj a 
-mv H — mu x H 

2 2 uj 



xvu — (v 2 + u 2 x 2 ) arctan \-^— 

\wx 



(5) 



The Lagrangian for the system (1) can now be constructed from the known 
expression [2] 

i( ,, l0 _f^i. (6 ) 
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For the general case (4), it is no possible to get a close expression for the 
Lagrangian, but for the weak dissipation case (5), one gets 



T 1 2 1 2 2, mw a 

L = -mv men x H 

2 2 lu 



(uj 2 x 2 — v 2 ) arctan ( — ) +ujxv In 

\UJXj 



' UJ 2 X 2 + v 2 ' 
uj 2 x 2 



The generalized linear momentum (p = dL/dv) is 
mui a 



p = mv + 



CO 



ujx + xuj In 



' ijJ 2 X 2 + v T 



oj 2 x 2 



— 2v arctan ( — ) 

UJX J 



(7) 
(8) 



To get the Hamiltonian, it is necessary from (8) to express the variable v 
as a function of the variables x and p, v = v(x,p). In this way, one makes 
the substitution of this variable on the Legendre transformation, H(x,p) = 
v(x,p) p — L(x, v(x,p)), or in the constant of motion, H(x,p) = K(x, v(x,p)). 
However, one notices immediately form (8) that it is not possible to do this. 
Thus, the Hamiltonian can not be given explicitly but implicitly through the 
constant of motion (5). 

A trajectory in the phase space (x, v) can be seen on Fig. 1, where the constant 
of motion (5) has been used. Of course, in order to keep the continuity at 
(v — 0,x < 0), the arctan function changes its value due to the multivalue 
functions. This means that our constant of motion is really a local constant 
of motion (which is valid on the half plane v < or v > 0), and it changes its 
value every time the trajectory crosses the line v — 0. This number of crossing 
is a numerable set. Therefore, this set has measure zero [10]. In this way, one 
can say that (4) or (5) represents a constant of motion almost everywhere in 
the phase space. 

Let us now change the variables (x, v) by a new variables (0, J) defined as 



<p = arctan — , 

yOJX, 



and J = — [v + 

2uj v 

where the inverse transformation is given by 

and 



2 2 

us x 



2J 



x = 



COS0 , 



2uJ . 1 
v — \ I sin . 



V muj V tu 

The dynamical system (2) is then transformed to the system 

2uj n tan < 



= -\uj + 



1 + tan 2 



and 



J = — 4u a J sin 2 . 
These equations are readily solved, and their solutions are 



<p(t) = — arctan 



t&nhvt + a) 

to 



(9a) 
(96) 

(10a) 
(106) 

(11a) 
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and 



J(t) = Jo e 



-u a f(t) 



(lib) 



where a and J Q are constants determinate by the initial conditions, and f(t) 
si defined as 



/(*) 



4(tanu;t — uj a /uo) t 



(11c) 



1 + (tanwt — uj a /uj) 

To know the trajectory in the new phase space (0, J), one can express J as a 
function of through the integration of J = (dJ/d(f))<fi. This brings about the 
following expression 



j(0) = J ° e ^arctan(tan0 + u a /uj) ? 
uj + u! a sin 



(12) 



where J Q is another constant. Fig. 2 shows Trajectories in this phase space for 
u) a = and for u a ^ 0, where the expected jumps at n/2 and 3n/2 are clearly 
seen at the scale shown. 



To finish the classical analysis, let us write the constant of motion, the La- 
grangian and the generalized linear momentum for the weak dissipative case 
in terms of the variables and J. These are given by 



K = ojJ + u a J 



sin(20) - 20 



L = — cjJcos(20) + 2uj a J 



cos(20) — 4 sin(20) ln(cos 0) 



and 



p = \JlmojJ sin + 



1 2m J 



cos 0(1 — 2 ln(cos 0) ) — 20 sin 



(13a) 



(136) 



(13c) 



Once again, one see the impossibility to get the Hamiltonian H(<f), J) due to 
complexity of the expression (13c). 



3. Quantization of the constant of motion 

Due to the impossibility of getting the Hamiltonian explicitly for the au- 
tonomous dynamical system (2), one may propose to extend the Shrodinger 
quantization to a dissipative system through the quantization of the constant 
of motion associated to it. This can be made by associating an Hermitian 
operator to the velocity as 

~ - (-\A\ 

m dx 

In this way, if K is the Hermitian operator associated to the constant of motion 
(which must have units of energy), the associated Shrodinger equation of the 
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classical autonomous system would be 

ih^ = K(x,v)V, (15) 

where \I/ = ty(x,t). Of course, the whole quantum mechanics structure is ex- 
actly the same but with the velocity operator instead of the linear momentum 
operator as the main operator of the quantum system. Since (15) represents 
an stationary problem, one just has to solve the eigenvalue problem 

K(x,v) <j) = E<j> , (16) 

where = <f)(x) and one has used ^(x, t) = 0(x) exp(—iEt/h) in (15). If K can 
be written as K = K Q + K I: where the solution of the problem K Q \n) = E^\n) 
is known, the eigenvalues of (16) are given at first order in perturbation theory 
[11] by 

E n = £f + (n\Kj\n) , (17) 

where one has used Dirac notation [12]. Our constant of motion (5) or 
(13a) can be expressed of the form K = K Q + Kj, where K Q = uoJ and 
Kj = uj a J(sin(20) — 20). It is well known that the eigenvalues of the har- 
monic oscillator without dissipation are given by E n °) = hu(n + 1/2). Thus, 
K Q is diagonal in the basis {\n)} and can be expressed in terms of ascent , a + , 
and descent, a, operators as 

K Q = Tiu(a + a + 1/2) , (18) 
with a and a + defined in terms of x and v as 



— - {lux + iv) , and a + = J — — (oux — iv) . (19) 
2noj V 2ncu 

These operator have the following commutation relations 

[a, a + ] = 1 , [a, a] = [a + , a + ] = , [a, a + a] = a , [a + , a + a] = —a + . 

Additionally, iV = a + a is called the number operator and is diagonal in the 
basis {|n)}, N\n) = n\n), This, in turns, implies that the operator associated 
to the variable J is given by 

J = h(N+ 1/2). (20) 

Thus, our main problem is to assign a Hermitian operator to the function 
Kj = 2cj a J(cos0sin0 — 0). According to refence [13], one has the following 
assignments 

sin — > S , cos — >■ C , — >■ , (21) 
where S and C are Hermitian operators such that 

S + uj 2 S = C + u 2 C = , [C, N]=iS , [S, N] = -iC (22a) 
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and 



§\ n ) = i(|„ + l>-| n -l>) ,C\n) = \ (W 1) + \n- 1>) ,[C,S} = ^, (22b) 

with 7r is the projector on the ground state of the quantum harmonic oscillator, 
7T D = |0)(0|. Therefore, C and S commute for any exited state. The operators 
4> is defined as 



= 



E 

k=0 



2k + l\ k J 



"1/2 \ ^2fc+l 



(23) 



where 



si the combinatorial factor. In this way, the operator associated 



to the function Kj can be written as 



JCS + CJS + CSJ + JSC + SJC + SCJ 



J(j> + (j>J 



• (24) 



The first term has not contribution at first order in perturbation theory be- 
cause of (22b), the action of the operators C or S increases or decreases the 
state number such that the expected value (n\..\n) is always zero. Due to the 
same reason, the contribution of the second term of Eq. (24) will come only 
from the term tt/2 of (23). Thus, at first order in perturbation theory, one has 
the following correction of the eigenvalues 

SEW = -hu a iv(n + 1/2) . (25a) 

Thus, one gets 

UJ a TT' 



E n w £<°> + 6E& = hu(n + 1/2) 



1 - 



That is, there is a small shift on the frequency of oscillation given by 



00 = U! 



■KUJo 



(26a) 



(266) 



At second order in perturbation theory,^ = E^+SE^+SE^\ the correction 
on the energy will come from the expression 



t'l" 1 771(0) 



Using (22b), (23) and (24), one has 



(k\Kr\n) = 



12 



(2k + 4n + 5)<J fc , n+2 - (2k + An + l)<5 fe , n _ 2 
21+1 (-1) 1 f-l/2\ (21 + 1 



(28) 



Therefore, the correction at second order can be written as 



6E?> 



,2 L 
—n H — 

3 4' 



00 21+1 

EE 

1=0 s=0 



-1/2 \ 2 [21 + l x 



(2n -21- s) 



(2l + lf(2l + l + s) 
(28) 
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Conclusions 



We have used the constant of motion for the one- dimensional dissipative har- 
monic oscillator to study the classical trajectories in the phase space and to 
point out the difficulty to get its Hamiltonian explicitly. Due to this problem, 
we have proposed the quantization of the constant of motion directly , via the 
association of the velocity operator. Then, the eigenvalues where calculated 
up to second order within perturbation theory to see the first effect of the dis- 
sipation on them. This effect at first order corresponds to have a shift by the 
quantity 7iu a on the frequency of oscillation of the nondissipative harmonic 
oscillator. 
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Figure captions 



Fig. 1 Trajectory on the phase space (x,v) as determinate by (4), where 
m = 1 Kgr, uj — 1 sec" 1 and u a = 0.001 sec" 1 . 

Fig. 2 Trajectory on the phase space ((f), J) as determinate by (12). The 
straight horizontal line corresponds^ to uj a = 0, solid line corresponds to 
u a = 0.001 sec" 1 , uj = 1 sec" 1 and J G = 1 Joules-sec. 
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